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Abstract 

We consider the evaporation of (4 + n) -dimensional non-rotating black holes into gravi- 
tons. We calculate the energy emission rate for gravitons in the bulk obtaining analytical 
solutions of the master equation satisfied by all three types (S,V,T) of gravitational 
perturbations. Our results, valid in the low-energy regime, show a vector radiation dom- 
inance for every value of n, while the relative magnitude of the energy emission rate of 
the subdominant scalar and tensor radiation depends on n. The low-energy emission rate 
in the bulk for gravitons is well below that for a scalar field, due to the absence of the 
dominant I = 0, 1 modes from the gravitational spectrum. Higher partial waves though 
may modify this behaviour at higher energies. The calculated low-energy emission rate, 
for all types of degrees of freedom decreases with n, although the full energy emission 
rate, integrated over all frequencies, is expected to increase with n, as in the previously 
studied case of a bulk scalar field. 



1 Introduction 



The quest for unification of fundamental interactions with gravity has led to the idea of 
extra spatial dimensions. String Theory, or some form of it, considered at present as a 
promising consistent theoretical framework for such a unification at the quantum level, 
requires extra spatial dimensions. In this higher-dimensional context, which has received 
a lot of attention in the last few years PQ |2] , gravity propagates in D = 4 + n dimen- 
sions (Bulk), while all Standard Model degrees of freedom are assumed to be confined on 
a four- dimensional D-Brane. In theories with large extra dimensions pQ the traditional 
Planck scale Mpi ~ 10 18 GeV is only an effective scale, related to the fundamental higher- 
dimensional gravity scale M* through the relation M Pl ~ M™ +2 R n , where R ~ (V n Y^ n is 
the effective size of the n extra spatial dimensions. If R ^> ip\ ~ 10~ 33 cm, the scale M* 
can be substantially lower than Mp\. If M* is sufficiently low, trans-planckian particle 
collisions could be feasible in accelerators or other environments. The product of such 
collisions would be objects that probe the extra dimensions. Black holes are among the 
possible products of trans-planckian collisions jSj- Black holes with a horizon size r# 
smaller than the size of extra dimensions R will be higher-dimensional objects centered 
at the brane and extending in the bulk. A classical treatment would require that the 
mass of the black hole Mpu would have to be larger than the fundamental gravity scale 
M*. Such black holes can be created in ground based colliders, although their appear- 
ance in cosmic rays is possible as well (for references, see the reviews jU E]). A black 
hole created in such trans-planckian collisions is expected to undergo a short "balding" 
phase in which it will shed its additional quantum numbers, apart from charge, angular 
momentum and mass. Then, after a more familiar Kerr-like phase, in which it will lose 
its angular momentum, a Schwarzschild phase follows during which the black hole will 
gradually lose its mass through the emission of Hawking radiation both in the bulk and 
in the brane, consisting of elementary particles of a very distinct thermal spectrum. This 
black hole radiation has been the subject of both analytical and numerical studies. This 
includes lower spin degrees of freedom |Hl El EE HD EE2j as well as graviton emission 
in the bulk [T3]. The study of graviton emission in the bulk corresponds to the study of 
perturbations in a gravitational background. The formalism for the treatment of gravita- 
tional perturbations of a higher-dimensional non-rotating black hole has been developed 
in jTl]. In the present article we consider the evaporation of (4 + n)-dimensional non- 
rotating black holes into gravitons and calculate the energy emission rate for gravitons in 
the bulk obtaining analytical solutions of the master equation satisfied by gravitational 
perturbations. Our results, valid in the low-energy regime, are complementary to exist- 
ing studies in the intermediate energy regime • These results show a vector radiation 
dominance for every value of n, while the relative magnitude of the energy emission rate 
of the subdominant scalar and tensor radiation depends on n. The low-energy emission 
rate in the bulk for gravitons is well below that for a scalar field, due to the absence of 
the dominant I = 0, 1 modes from the gravitational spectrum, although higher partial 
waves are likely to modify this behaviour at higher energies. The calculated low-energy 
emission rate, for all types of degrees of freedom decreases with n, although the full 
energy emission rate, integrated over all frequencies, is expected to increase with n, as 
in the previously studied case of a bulk scalar field. 
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2 Theoretical Framework 



According to the Large Extra Dimensions scenario PP, our universe is made up of the 
usual 3+1 non-compact dimensions plus n additional, compact, spacelike ones. The grav- 
itational background around a spherically-symmetric, neutral black hole formed in such 
a (4 + n)-dimensional, flat spacetime is described by the following line-element [13] 



ds 2 = -f(r)dt 2 + ^ r) +r 2 dQ 2 2+n , (1) 



with the metric function f(r) given by 



n+l 



f(r) = 1 - [f ) , (2) 

and dQ 2 +n denoting the area of the (2 + n)-dimensional unit sphere. The above black 
hole is characterized by a non-vanishing temperature given by 

^ n + 1 

Th = ~. , 3 

Attth 

and therefore emits Hawking radiation [6 with an energy emission rate that resembles the 
one for a black body. The strong gravitational field surrounding the black hole modifies 
the spectrum through the modification of the corresponding absorption probability \A\ 2 , 
that now depends on the energy u and spin s of the emitted particle as well as on the 
number of extra dimensions n (see, for example, [3]). The absorption probability can be 
found by solving the equation of motion of a particular field propagating in the vicinity of 
the black hole and by using classical scattering theory. As mentioned in the Introduction, 
in this work, we focus on the decay of a (4 + n) -dimensional, Schwarzschild black hole 
through the emission of gravitons in the bulk. 

Following a graviton in D dimensions can be decomposed into a symmetric 
traceless tensor, a vector and a scalar part. These are further expanded in terms of the 
spin-weighted spherical harmonics on the S n+2 unit sphere. The radial parts of all three 
types of gravitational perturbations are found to satisfy a Schrodinger-like equation of 
the form ^3] 

d ( „ d§\ , 2 



v - /(r 



M ; *J +(W ,/)4,=0 - (4) 



The potential V has a different form for each type of perturbation, namely 

n{n + 2) k(n + 2) 2 (r H \ n + l ~ 

r J 

for tensor- and vector-like perturbations, with k = —1 and k = 3 respectively, and 



;(i+n+1) + »'» + 2 )- <: (» + 2 ) 2 (; 



(5) 



f( r ) qx 3 + px 2 + wx + z 

KS — n 9 ) (6) 

r 2 4[2m+ (n + 2)(n + 3)x} 2 



for scalar gravitational perturbations, with 

n+l 



m = l(l + n + l)-n-2, x = J =1-/, (7) 
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and 



q = (n + 2) 4 (n + 3) 2 , z = 16m 3 + 4m 2 (n + 2)(n + 4) , 

p = (n + 2)(n + 3) [4m (2n 2 + 5n + 6) + n(n + 2)(n + 3)(n - 2)] , (8) 

w = -12m (n + 2) [m(n - 2) + n(n + 2)(n + 3)] . 

3 Solving the Field Equations 

In this section, we proceed to solve analytically the above equation for all three types of 
gravitational perturbations in the vicinity of a (4 + n)-dimensional neutral, spherically- 
symmetric black-hole. For this purpose, we will use a well-known approximate method 
and solve first the aforementioned equation at the two asymptotic radial regimes: close 
to the black hole horizon (r ~ r H ), and far away from it (r 3> r H ). The two solutions will 
then be stretched and matched in an intermediate zone to create a smooth analytical 
solution extending over the whole radial regime. 



3.1 The Near-Horizon Regime 

Here, we focus on the derivation of the solution for all types of gravitational perturbations 
in the asymptotic regime close to the black hole horizon. Due to the similar form of 
their potential, tensor and vector perturbations will be treated together, while the scalar 
perturbations will be dealt with separately. 

A. Tensor and Vector Perturbations 



We start by making the change of variable r — > f(r), in terms of which the field equation 
for tensor and vector perturbations can be written in the form 



d 2 $ 

/(W)^ + 



(2n + 3) 



./ 



(n + 1) 

(cor H ? 



+ 



k(n + 2) 



(n + l) 2 /(l-/) (I"/) 4(n+l) 



$ = 



(9) 



In the above, k — — 1 for tensor perturbations and k = 3 for vector ones, and we have 
also defined for convenience the quantity 



A _ 1(1 + n + 1) n(n + 2) 



(n + 1) 2 4(n+l) 2 



(10) 



If we now make the following field redefinition: $(/) = f a {l — f) l3 F(f) ) the above 
equation takes the form of a hypergeometric equation 



fa-f)^ + [c-(l + a + b)f]^-abF 



0, 



;n) 
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under the identifications 



(n 4- 2) (n 4- 2) 

* = u + P + ^T7S + G > b = a + P+^--^~G, c=l + 2a, (12) 
2(n + 1) 2(n + 1) 

where G an arbitrary constant. Demanding further that the coefficient of F(f) in our 
equation be indeed — ab yields three more additional constraints that determine the 
remaining unknown constants, a, /3 and G. Their corresponding values are found to be 

Q{T> v) = (l + fc)(n + 2) = ± iwrn_ 

4(n + l) n + 1 

and 

P = 2(^1) { _1 ± V (2Z + n + 1)2 " 4 ^} • (14) 

The indices of the hypergeometric equation (a, b, c) are now fully determined, which 
allows us to write the general solution of the hypergeometric equation, and thus the 
solution in the near-horizon regime, as 

$ NH (f)=A l r(l-ffF(a,b,c; f) 

+ A 2 r a (l-fyF(a-c+l,b-c+l,2-c;f), (15) 

where A\ 2 are arbitrary integration constants. There is, however, a boundary condition 
that this general solution must satisfy: no outgoing waves must be found right outside 
the black-hole horizon, as nothing can escape from within this area. To ensure this, we 
expand our solution in the limit r — > r#, or equivalently / — > 0, in which case we find 

$ NH (f) ~A 1 f a + A 2 f- a = Ai e- wy + A 2 , (16) 

where, in the last part of the above equation, we have used the 'tortoise-like' coordinate 

d JL = ^ . (17) 

The asymptotic solution is therefore written in terms of incoming and outgoing plane 
waves, as expected, since very close to the horizon the potential V for all types of grav- 
itational perturbations vanishes. Nevertheless, the aforementioned boundary condition 
forces us to discard the outgoing wave by setting A 2 = 0, that brings the near- horizon 
solution to its final form 

$N H (f) = A x r (1 - ff F(a, b, c; /) . (18) 



The only task remaining is to fix the arbitrary signs appearing in the expressions of 
both a and (3. As we may see from Eq. ([16)1 . the interchange a + <-» a_ would simply 
interchange the integration coefficients A\ <-> A 2 , therefore the sign of a can be chosen at 
random; here, we have chosen a = a_. On the other hand, the sign in f3 can be fixed by the 
convergence condition of the hypergeometric function, i.e. Re (c—a—b) = —-^^ — 2(3 > 0, 
that clearly demands that we choose f3 = /3_. 
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B. Scalar Perturbations 



By employing the same change of variable, the field equation for scalar gravitational 
perturbations can be brought to the form 



/(W)^ + 



(2n + 3) 



(n + 1) 

(ur H ) 2 



d<S> 



(n + 1) 2 /(1 - /) 16(n+ l) 2 m 2 (l - /) 



C 



$ = 0, (19) 



where we have defined the quantity 



C = 



q (l-ff+p(l-f)+w 
4(n + l) 2 [2m + (n + 2)(n + 3)(1 - f)f 



with 



p = p 



z(n + 2) 2 (n + 3f 
Am 2 



w 



w 



z(n + 2)(n + 3) 



m 



(20) 



(21) 



oo, 



In this form, Eq. (fT9*j) has poles at / = and / = 1 (or, at r = Th and r 
respectively) while the quantity C takes on a constant value in both limits. 

We then follow a similar method as before, and bring again Eq. (|T9*j) to the form of 
a hypergeometric equation, with the indices (a, b, c) given by Eqs. (fT2Jl. and the powers 
(a, j3) by Eqs. (|13|) - (jl4|) . The only difference arises in the value of the arbitrary constant 
G that, in the case of scalar perturbations, takes the value 



G (S) 



1 



2(n + l) 



(n + 2) 



q + p + w 



[2m + (n + 2)(n + 3)] i 



(22) 



By applying the same boundary condition of no-outgoing waves near the horizon, 
the general solution for scalar perturbations in the near-horizon regime is given again by 
Eq. (fTHj). with p = P_ and a = a_, as before. 



3.2 The Far-Field Regime 

We now turn to the far-field regime. In the limit r » r^, / -> 1, and the field equation for 
all types of gravitational perturbations - tensor, vector and scalar - takes the simplified 
form 

£ + (»>-^).-o. 

where A was defined in Eq. fllOjl . By further setting $ = y/r R, the new radial function 
R is found to satisfy a Bessel differential equation 

«"* + i« + f 1 _4 , U = 0, (24) 



dz 2 z dz \ z 2 
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with z = ojr and v — I + (n+ l)/2. Therefore, the general solution of Eq. ()23|). standing 
for the analytical solution in the far-field regime, can be written as 

$ffO) = Biy/r Jz+(„+i)/2 (wr) + 5 2 \/r Vj+(n+i)/2 (^) , (25) 

where J u and Y v are the Bessel functions of the first and second kind, respectively, and 
Bi t 2 arbitrary integration constants. 



3.3 Constructing the complete solution 

Having derived the two asymptotic analytical solutions near the black-hole horizon and 
infinity, we may now proceed to construct a solution valid at the whole radial regime by 
smoothly connecting the two solutions at an intermediate point. To this end, we need to 
stretch the near-horizon solution (JTSj) towards large values of r, and the far- field solution 
(|25|) towards small values of r. 

Before, however, being able to do the former, we need to rewrite the near-horizon 
solution in an alternative form where the argument of the hypergeometric function has 
been shifted from / to 1 — /. To this end, we use the standard formula [TH] 

F(a,b,c;f) = r(c) T(c - a - b) b +b _ c+11 _ f) (26) 
1 (c — a) 1 (c — b) 

+ (1 - fr a - b r(c) r r , ( ° + * ~ C) F(c - a, c - b, c - a - b + 1; 1 - /) 

r(a)r(6) 

in Eq. (|18|) . and subsequently take the limit r ^> rjj, or / — > 1. We then find 
x . , „ /r^\/3(«+i) r(c)r(c-a-fo) A / rjf \-i-0(n+i) r(c)r(a + 6-c) 

•™M - a, (_) r< /_^ )r(c _ t + 1 ( T ) r ( a)r W • 

(27) 

We now turn to the FF solution, which we expand in the opposite limit of r — ► 0. 
By using standard formulae for the Bessel functions ^H], we then obtain 



Although we have brought both asymptotic solutions in a power-like form, the powers 
are not the same. One way to simplify the matching procedure is to take the low-energy 
limit u)rji <C 1 in the expression of (3, Eq. (JHJ). At first-order approximation, the ujth- 
term may be ignored, in which case the powers in Eqs. (|2"7j) and (f^Hj) become identical. A 
smooth matching is then achieved, and a complete solution is constructed, if the following 
relations between the near-horizon and far-field integration constants hold 



Bi _ ( 2 Y+— r(l + ^) T(c)T(c-a-b) 
M \ur H ) r(c - a) T(c - b) 



(29) 



(^H\ l+ ^ T(c)T(a + b-c) 

M 71 \ 2 J r(/ + H±i) r(a)r(6)Vri" 



6 



The above completes the derivation, in an analytical way, of the solution for all types 
of gravitational perturbations in a (4 + n)-dimensional Schwarzschild black hole back- 
ground in the low-energy regime. We now proceed to the calculation of the corresponding 
absorption coefficient. 



4 The Absorption Coefficient 

The far-field solution (J25|) can also be expanded in the limit r — * oo, where the effective 
potential for all types of gravitational perturbations again vanishes due to the asymp- 
totically flat behaviour of the metric. In this asymptotic regime, we therefore expect the 
general solution to be described again by incoming and outgoing plane waves. Indeed, in 
the limit r — > oo, we obtain 

$ FF (r) ~ -^L= { (B 1 - iB 2 ) e<^-5 + (B x + iB 2 ) ) + ..., (31) 

\i2txuj *> J 

where as before v — I + (n + l)/2. From the above expression, we can easily define the 
reflection coefficient TZi as the ratio of the amplitude of the outgoing wave over the one 
of the incoming wave. Then, the absorption probability follows from the relation 

\Ai\ 2 = l- |7^| 2 = 1- 

where B is defined as 

b 1 = _ f j_\ 2l+n+1 r(nf)r(i + f) r(q)r(6)r(c-a-6) 

B 2 \tor H J 7rT(c-a)T(c-b)T(a + b-c) ' 1 ' 

The above two equations comprise our main analytical result for the absorption co- 
efficient associated with the propagation of gravitons in the higher- dimensional black 
hole background given in Eq. Individual solutions for scalar, vector and tensor type 
of perturbations may easily follow upon substituting the corresponding values for the 
hypergeometric indices (a, b, c) found in the previous section. 



B 



B + i 



(32) 



4.1 Simplified analytical result 

The expression can be simplified, and a compact analytical result for the absorption 
coefficient may thus be derived, if we further expand it in the low-energy limit urn 1. 
For convenience, we re-write the hypergeometric indices as 

a = a + (3 + d, b = a + (3 + G 2 , c=l + 2a, (34) 

where ^ ^ 

Gl ^2l^Ti) +G ' G2 ^WTT)~ G - (35) 
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As G takes a different value for scalar, vector and tensor type of gravitational perturba- 
tions, G12 will also depend on the type of perturbation studied. We also write Eq. (J32J) 
in the form 



2i(B* - B) 

1^1 " BB* + i(B* -B) + 1' {6b) 

and note that, in the low-energy limit, BB* ^> i(B* — B) ^> 1. Therefore, by keeping 
only the dominant term in the denominator, we arrive at 



\A i \ 2 = K{ur H ,P) {Z(a,f3)-Z*(a,f3)} , (37) 



where 



irt a\ (ur„\*+*M 2m (I + 2±1) r(-l + 2B + G 1 + G 2 ) 

*(«*■„/» s - (— ) r(t + ! ^y r(i-2g-G,-G,) (38) 

and 

7 , m _ r(l +a-(3- G^rq + q-^-Ga) 

Let us focus first on the expression for K(ur H , [3). By using the following identity 
satisfied by the Gamma functions [21] 

T{x)T(l - x) = -xT(-x)T(x) = — ^ , (40) 

sin(7rx) 

keeping the dominant term (3^ in the expansion of (3 (Eq. ITU) in the limit uth -C 1, i.e. 



(2/ + , + 2) 

P 2(n+l) V ; 



and using (jHHJ), we find 



/cjr // \ 2/+n+1 m 2 (n+l) 
K(urH, p) = —p— — -tto - x2 , ; 



V 2 / 



r (Z + r (1 + sin [7r(2/5<°> + d + G 2 )] 



(42) 



The only complex quantity appearing in the expression of Z(ot,/3) is a, which is 
purely imaginary, therefore Z*(a,f3) = Z(—a,f3). By using again Eq. flffij), we may write 

7T 2 

^ ~ Z * = \T(a + (3 + G 1 W\T(a + (3 + G 2 )\ 2 (43) 

sin [7r(a + (3 + Gi)] sin [7r(a + (3 + G 2 )] — sin [7r(a — (3 — Gi)] sin [7r(a — /3 — G 2 )] 
X | sin [7r(a + (3 + G x )] | 2 | sin [vr(a + /? + G 2 )] | 2 ' 

Expanding again in the limit ujth <C 1, or equivalently a — > 0, we obtain 

Z - Z* = — sin \u(2(3^ + G x + G 2 )l T(l - /3 (0) - G x ) 2 T(l - /3 (0) - G 2 ) 2 . (44) 

7T 



s 



a(T)\2 |/l( V )|2 |/l( S )|2 

*/t.j I I j I I 



n 


= 2 


/ = 2 


1.7 


• 10 


-4 




7.4 


• 10 


-3 


(ur H ) 8 


2.7 


• 10- 3 (ujr H ) s 






Z = 3 


1.1 


10" 


-6 


(urn) 10 


1.6 


io- 


- 5 


[urn) 10 


1.8 


10- 5 (cur H ) 10 






I = 4 


4.6 


io- 


" 9 


i^r H ) 12 


3.7 


io- 


8 




5.5 


10- 8 (ujr H ) 12 


n 


= 4 


/ = 2 


3.2 


io- 


- 6 


(uJr H ) 10 


2.2 


io- 


1 


[uJr H ) 10 


3.4 


10- 5 (ujr H ) w 






Z = 3 


1.9 


io- 


-8 


(u;r H ) 12 


4.2 


io- 


- 7 




3.1 


lQ- 7 {ur H ) 12 






Z = 4 


8.1 • 


10" 


11 


(ujr H ) u 


9.5 • 


10" 


-10 


(ujr H ) u 


1.3 


■ 10- 9 (tur H ) 14 


n 


= 6 


Z = 2 


3.1 ■ 


10" 


8 


'ur H ) 12 


3.2 


10" 


6 


[uJr H r 


2.2 


10- 7 (ior H ) 12 






Z = 3 


1.5- 


io- 


10 




4.7 


io- 


- 9 




1.9 


10- 9 {ur H ) M 






Z = 4 


5.4- 


io- 


13 


K) 16 


8.9- 


10" 


-12 


{ur H ) 16 


8.4- 


10- 12 (iur H ) w 



Table 1: Dependence of the absorption probability for tensor, vector and scalar gravita- 
tional perturbations in the bulk on n and Z, in the asymptotic regime ujth — > 0. 



By putting Eqs. (J42j) and (|44|) together, we finally obtain the following expression 
for the absorption probability in the asymptotic low-energy regime 



\M 2 



4tt ( 



ur H Y 



r(z + n±3) 2 r(i + ^)' 



(45) 



In the above, we have used again the zero-order approximation ()41)1 for /3, and the 
definitions (jHSJ) to recover the dependence on the parameter G. According to Eq. (JTHJ), 
the value of the latter parameter is zero for tensor gravitational perturbations and ^| 
for vector ones, while for scalar gravitational perturbations its value is given in Eq. 
(|2*2*j). In the case of tensor perturbations, it can be shown that the above result reduces 
to the one for the absorption probability for a scalar field propagating in the bulk [7j. 
This result should have been anticipated by looking at the equation satisfied by the 
tensor-like gravitational perturbations in the bulk: starting from Eq. (J1J) and setting 
$Cr)(r ' 



V r n+2 $(r), the new radial function $(r) is found to satisfy the equation 



fir) d 

r n+2 ^ r 



dr 



+ 



^Z(Z + n + l) 



$(r) = 



(46) 



that describes indeed the propagation of a scalar field in the black-hole background of 
Eq. © 0. 

From Eq. ([45)1 . one may see that the absorption probability depends on both the 
angular momentum number Z and number of extra dimensions n, through the arguments 
of the Gamma functions as well as the power of ujth- As either Z or n increases, the 
latter increases too, which, for ur H -C 1, causes a suppression in the value of \Ai\ 2 . As 
the behaviour of the remaining factor is not equally clear, in Table 1 we display the 
explicit value of \A{\ 2 for all three types of gravitational perturbations, as these follow 
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ujth \A\ 1 2 (simplified expression) \A\ | 2 (exact expression) 

0.001 1.6997 x 10" 28 1.6997 x 10" 28 

0.01 1.6997 x 10~ 20 1.6999 x lO" 20 

0.1 1.6697 xlO^ 12 1.7112 x 10~ 12 

0.3 1.1152 x 10~ 8 1.1839 x 10~ 8 

0.5 6.6396 x 10~ 7 7.8001 x 10" 7 



Table 2: Deviation between the values of the absorption probability given by the simpli- 
fied and complete analytical expression, for tensor gravitational perturbations, for n = 2, 
I = 2 and different values of uor^ < 1. 

from the simplified expression (J4*5*j) . for the indicative values n = 2,4, 6 and I = 2,3,4. 
From these entries, one may easily conclude that, as either I or n increases, the value of 
the absorption probability for all types of gravitational perturbations in the asymptotic 
low-energy regime is significantly suppressed. 

As a final comment, let us note here that, for the same values of / and n, \A{\ 2 assumes 
a different value for different types of gravitational perturbations. From the entries of 
Table 1, we may easily see that the tensor perturbations are by orders of magnitude 
suppressed compared to the vector and scalar perturbations, while the relative magnitude 
of the latter two is strongly dependent on the particular values of I and n. As the tensor 
perturbations have the same absorption probability as the bulk scalar fields, it would be 
interesting to see whether gravitons dominate over scalar fields during the emission of 
Hawking radiation in the bulk. 

4.2 Complete analytical result 

In deriving our main result for the absorption probability for gravitational perturbations 
in the bulk, Eqs. (|3^ j) - (j3^j) . the use of the low-energy assumption was made only once 
- during the matching of the two asymptotic solutions in the intermediate zone. Nev- 
ertheless, that was enough to restrict the validity of our solution to values of the uoth 
parameter well below unity. The simplified analytical result (|43|) was, on the other hand, 
the result of a series of expansions in the arguments of the Gamma functions appearing 
in Eq. (}3*3*j) . and, as a result, its validity is significantly more restricted. In Table 2, we 
display the values of the absorption coefficient derived by using our two analytical ex- 
pressions, Eqs. fl3^|) - (jHH|) and (j43|) . as ur H ranges from 0.001 to 0.5. We may easily see 
that, for very low values of urn, the agreement between the two values is remarkable; 
however, as soon as urn reaches the value 0.5, the deviation between the two values 
appearing in the last row of Table 2 reaches the magnitude of 15%. 

In this section, we return to our main analytical result for the absorption probability 
of gravitons in the bulk, Eqs. (}3*2*|) - (j3*3*j) . As it was shown in the case of emission on the 
brane [3JE], the analytical, non-simplified, expression for \A{\ 2 is in excellent agreement 
with the exact numerical result in the low-energy regime and in a good - both quantita- 
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Figure 1: Absorption probability \Ai\ 2 for tensor (solid lines), vector (short-dashed lines) 
and scalar (long-dashed lines) gravitational perturbations in the bulk for: (a) n — 6, and 
I = 2, 3, 4, 5, and (b) I = 2, and n = 0, 2, 4, 6. 



tive and qualitative - agreement in the intermediate-energy regime. In the high-energy 
regime, the validity of our expression naturally breaks down. There, we expect the grey- 
body factor (absorption cross-section) to be given in terms of the absorptive area of the 
black hole: this result has been computed in EI , and holds for all types of perturbations 
independently of their spin. 

In Figs. 1(a) and (b), we depict the absorption probability for all types of gravita- 
tional perturbations, as a function of the dimensionless energy parameter uth, and for 
different values of the angular momentum number I and number of extra dimensions 
n, respectively. As expected, for all types of perturbations and values of I and n, the 
absorption coefficient vanishes when the energy of the propagating particle goes to zero, 
while it increases with ujth- This is in agreement with the classical scattering theory: 
the larger its energy, the more likely it is for a particle to overcome the gravitational 
field of the black hole and escape to infinity. Figure 1(a) reveals that, as I increases, 
the absorption probability for all types of gravitational perturbations is suppressed, in 
accordance with the behaviour noticed in the entries of Table 1; although the results 
depicted correspond to the case n = 6, this behaviour was found to hold for all values 
of n. For fixed on the other hand, as we may see from Fig. 1(b), \Ai\ 2 is found to 
exhibit a significant suppression also with n, for all types of gravitational perturbations, 
in agreement again with the behaviour noticed in the entries of Table 1. This behaviour 
is observed for the lowest two partial modes with I = 2 and / = 3; for higher values of 
I, \Ai | 2 shows a temporary enhancement, as n increases from to 2, that nevertheless 
changes again to a rapid suppression as n increases further. 

As it was noticed by studying the asymptotic low-energy regime, the absorption 
probability for tensor gravitational perturbations is significantly suppressed, compared 
to the other two types, also in the extended low- and intermediate-energy regime. On 
the other hand, the vector and scalar perturbations may dominate one over the other 
for different values of /: as it can be seen from our figures, the vector-like gravitational 
perturbations predominantly dominate for low values of /, while the scalar-like perturba- 
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48 


2772 


25 


450 


I 


= 5 


64 


23868 


70 


7140 


36 


1122 



Table 3: Multiplicities of states corresponding to the same angular momentum number I 
for tensor, vector and scalar perturbations, for n = 1 and n = 6. 

tions usually take over for large values of /. As the higher partial waves are significantly 
suppressed, we are led to conclude that the vector-type perturbations will be the dom- 
inant gravitational degree of freedom to be emitted in the bulk by a Schwarzschild-like 
higher- dimensional black hole. However, different types of perturbations are character- 
ized by a different multiplicity of states for the same value of /, and the latter must be 
taken into account before such conclusions can be safely drawn. 



5 Energy Emission Rate 

We finally turn to the energy emission rate by the higher-dimensional Schwarzschild 
black hole for gravitons in the bulk. Let us first address the issue of the multiplicities 
of states that correspond to the same angular momentum number I. For tensor, vector, 
and scalar type of perturbations these are given by the expressions [17] 

(T) _ njn + 3)(l + n + 2)(l-l)(2l + n + l)(l + n- 1)! 

1 2(Z+l)!(n+l)l ' 1 } 

{v) l(l + n + l)(2l + n + l)(l + n-l)\ 

Nl = (TTT)!^! ' (48) 

{s) _ (21 + n+m + r^ 

1 " l\(n+l)\ ' 1 } 

respectively. As it is evident, the multiplicities depend also on the dimensionality of 
spacetime. In Table 3, we display the multiplicities of states for tensor, vector and scalar 
perturbations, for some indicative values of / and n. We immediately observe the prolif- 
eration of states as either one of these two parameters increases, especially the one with 
the number of extra dimensions. It thus becomes clear that the value of the absorption 
probability alone will not be the only decisive factor that will determine the contribution 
of each type of gravitational perturbation to the total emission rate of the black hole. 

With the multiplicities of states and values of the absorption probability for each 
type of perturbation at our disposal, we can now calculate the corresponding energy 
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Figure 2: Energy emission rates for tensor (thin solid lines), vector (short-dashed lines) 
and scalar (long-dashed lines) gravitational perturbations, and scalar fields (thick solid 
lines) in the bulk, for (a) n — 1, and (b) n — 6. 



emission rate at the low-energy regime. The contribution of each type of gravitational 
perturbation to the total graviton energy emission rate is given by the expression [3] 



where the superscript P = (T, V, S) denotes the type of perturbation, and the black hole 
temperature is given in Eq. (|3*|). Summing over the three contributions, we may obtain 
the total amount of energy emitted per unit time and unit frequency by the black hole in 
the form of gravitons in the bulk. As this result would be more useful in the context of 
an exact numerical analysis (to which we hope to return soon), here we will concentrate 
on the relative emission rates for the different types of gravitational perturbations and 
their relative magnitude to the one for bulk scalar fields. 

A simple numerical analysis, by combining Eq. ([50)1 with the entries of Tables 1 and 
3, reveals that, in the asymptotic low-energy regime, the vector-like perturbations are 
indeed the dominant type of gravitational degree of freedom emitted in the bulk by the 
black hole. For example, for n = 2 and I = 2, vector perturbations amount to 85% 
of the total gravitational degrees of freedom emitted, compared to 13% for scalar and 
2% for tensor degrees. As n increases further, so does the dominance of the vector-like 
perturbations that, for n = 6 and I = 2, reaches the magnitude of 97%. This dominance 
is never over-turned but it may be significantly decreased, at the level of higher partial 
waves: for instance, for n = 6 and I = 4, the vector, scalar and tensor perturbations 
correspond to 74%, 12% and 14%, respectively, of the total number of gravitational 
degrees of freedom emitted. In fact, due to their large multiplicity of states, the tensor 
perturbations dominate over the scalar ones, for most large values of n and/or I. 

The above results valid at the lowest part of the radiation spectrum may change 
for higher values of the energy parameter. Therefore, in Figs. 2(a) and (b), we plot the 
energy emission rates for different types of gravitational perturbations, for n = 1 and 
n = 6, respectively, as these follow by using our complete analytical expression (p ^ - lp^j) 




(50) 
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for the absorption probability. For comparison, we also display the energy emission rate 
for scalar fields in the bulk. In the sum over I in Eq. (J5Uj) we have included all modes up 
to I = 12, although at the very low-energy part of the spectrum, the contribution of all 
modes with / > 4 is at least four orders of magnitude smaller than the one of the / = 2 
for every value of n. Having included a large enough number of modes in the sum, allows 
for deviations from this behaviour as the energy increases. 

The results depicted in Fig. 2 are in fact in excellent agreement with the conclusions 
derived by using the simplified expression for \Ai\ 2 . The vector- type perturbations are 
indeed the dominant gravitational degrees of freedom to be emitted by the black hole 
in the bulk for every value of n. What depends on the number of extra dimensions is 
the relative magnitude of the energy emission rate of scalar and tensor perturbations: 
whereas for low n, the tensor perturbations are subdominant to the scalar ones over 
the entire low-energy regime, for high n, they clearly dominate over the latter ones. 
Despite the above, the energy emission rates for all types of gravitational perturbations 
in the bulk, even when combined, remain well below the one for scalar fields in the low- 
energy regime. This is due to the fact that the emission rate for scalar fields receives 
a significant enhancement, in this particular energy regime, from the dominant I = 
and I = 1 modes, that are absent from the spectrum of gravitational perturbations. As 
the energy parameter increases though, we expect the higher partial waves to gradually 
come into dominance and possibly help gravitons to dominate over the bulk scalar fields. 

By comparing finally the vertical axes of the two plots in Fig. 2, we conclude that 
the low-energy emission rate for all types of degrees of freedom in the bulk decreases 
as the number of extra dimensions increases. Although both the temperature of the 
black hole and multiplicity of states undergo a significant enhancement as n increases, 
the equally significant suppression of the absorption probability, depicted in Fig. 1(b), 
prevails, leading to the suppression of the number of degrees of freedom emitted by 
the black hole. This low-energy suppression for bulk scalar fields was first seen in [7J |5] 
but the exact numerical analysis performed in the latter work showed that, for higher 
values of the energy parameter, the spectrum is enhanced with the number of extra 
dimensions. This is caused by the milder suppression of the absorption probability with 
n at higher frequencies, but also by the shift of the emission curve towards higher energies, 
in accordance with Wien's law - the latter leads to the emission of more higher frequency 
particles and less low-energy ones as n increases [9 . Due to the similarities observed in the 
behaviour of gravitational and scalar fields in the bulk, we expect the same enhancement 
to take place also for gravitons at higher energy regimes. 

6 Conclusions 

A higher-dimensional black hole created in the context of a brane-world theory will 
decay through the emission of Hawking radiation both in the bulk and on the brane. 
The type of particles emitted along each "channel" are determined by the assumptions 
of the particular model, with only gravitons and possibly scalars propagating in the bulk 
within the framework of the Large Extra Dimensions scenario. Although the emission of 
scalars has been studied in detail, until recently gravitons have received little attention. 
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In this work, we have addressed this gap in the literature by investigating the emission 
of tensor, vector and scalar gravitational modes in the bulk from a (4 + n)-dimensional 
Schwarzschild black hole. Working in the low energy regime, we analytically solved the 
corresponding field equations and computed the absorption probability in each case. Both 
a complete analytic expression and its asymptotic low-energy simplification were studied 
in detail, and their dependence on the angular momentum number / and number of extra 
dimensions n was examined. Although numerically different as the energy increases, these 
two expressions have identical qualitative behaviours, that reveal an increase in the value 
of the absorption probability with increasing energy and suppression as the number of 
partial wave I and extra dimensions n increase. 

The complete analytical expression for the absorption probability was then used to 
derive the contribution of each gravitational degree of freedom to the total graviton 
emission rate of the black hole in the bulk. Accounting for the rapid proliferation of 
state multiplicities with I or n, a sum over the first twelve partial waves was performed 
to obtain the total low-energy emission rate for each gravitational degree of freedom. Our 
results show that vector perturbations are the dominant mode emitted in the bulk for 
all values of n. The relative emission rates of the subdominant scalar and tensor modes 
depend on n, with scalars foremost at small n and tensors more prevalent at high n. 
The absence of the I = 0, 1 partial waves, dominant in the low-energy regime, from all 
gravitational spectra causes even the total gravitational emission rate to be subdominant 
to that of scalar fields in the bulk. Finally, as previously found for bulk scalar fields, the 
energy emission rates for all types of gravitational perturbations are suppressed with the 
number of extra dimensions in the entire low-energy regime. 

In this work we have focused on the low-energy part of the Hawking radiation spec- 
trum as it permitted analytical calculation and derivation of closed form expressions for 
the absorption probability for gravitons. Although the radiation emitted in the bulk is 
not directly observable, it determines the energy left for emission on our brane. In this 
context, our results, in addition to their theoretical interest, would be of particular use 
to experiments developed to detect the low-energy spectrum of radiation emitted from 
a higher- dimensional black hole. As the exact form of the complete gravitational spec- 
trum is still pending, we hope to return in the near future with results from a numerical 
analysis that could only provide the answer to this question. 

Note added in proof. While this work was at its last stages, two relevant papers 
appeared in the literature, JH] and PH] • in particular, the latter work [T.9., , that studies the 
graviton emission rate in the bulk, overlaps with ours with respect to the analytical results 
for the cases of tensor and vector gravitational perturbations, which are in agreement 
with ours. 
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